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. . . We investigate how the dynamical production of quantum entanglement for weakly coupled, com- 

CS| ■ posite quantum systems is influenced by the chaotic dynamics of the corresponding classical system, 

' using coupled kicked tops. The linear entropy for the subsystem (a kicked top) is employed as a 

, measure of entanglement. A perturbative formula for the entanglement production rate is derived. 

The formula contains a correlation function that can be evaluated only from the Information of 
^ ' uncoupled tops. Using this expression and the assumption that the correlation function decays ex- 

Q , ponentially which is plausible for chaotic tops, it is shown that the increment of the strength of chaos 

■ does not enhance the production rate of entanglement when the coupling is weak enough and the 

' subsystems (kicked tops) are strongly chaotic. The result is confirmed by numerical experiments. 

^\ I The perturbative approach is also applied to a weakly chaotic region, where tori and chaotic sea 

coexist in the corresponding classical phase space, to reexamine a recent numerical study that sug- 
gests an intimate relationship between the linear stability of the corresponding classical trajectory 
and the entanglement production rate. 

> ' 

O ■ PACS numbers: 05.45.Mt,03.65.Ud,05.70.Ln,03.67.-a 



CN ■ I. INTRODUCTION 

o . 

Quantum entanglement (in short, entanglement) in composite systems has been discussed as a paradoxical issue [Q, 
but it is experimentally confirmed, and utilized in quantum Information processing Although the definition of 
entanglement itself is not of dynamical nature, entangled states are often generated dynamically ^, ^. That is, 
even if subsystems are not entangled initially, the interaction between them produces entanglement in the system 
as time elapses. It is easily expected that the dynamical production of quantum entanglement heavily depends on 
the qualitative nature of dynamics. An important qualitative distinction of quantum dynamics is provided by the 
IL* ' corresponding classical dynamics, whether it is regular or chaotic, as is well known in the literature of "quantum 
.5^ ' chaos" 1^. Hence, there have been investigations to answer the foUowing question |^, ^, ^, Does the production 
of entanglement depend on whether the corresponding classical system is chaotic or regular? The authors of Refs. j^, 
5-H 01 Hi ü concluded that chaotic swtems tend to produce larger entanglement than the regular systems in general 
^ (exceptions are shown in Refs. |0, |lï|). Accordingly, it is natural to raise the next question: For chaotic systems, 
does the strength of chaos increase the degree of entanglement? This issue was first investigated by Miller and 
Sarkar jï^. They employed coupled kicked tops (CKTs) as their model system, and numerically found that the von 
Neumann entropy of the subsystem linearly depends on the sum of positive (finite-time) Lyapunov exponents of the 
corresponding classical systems. To authors' knowledge, however, there has been no theoretical explanation for their 
result. 

In this paper we examine the same system (CKTs) to elucidate the mechanism of the entanglement production. In 
particular, we clarify how dynamical aspects of quantum entanglement for CKTs are affected by the chaotic dynamics 
of the corresponding classical CKTs when the coupling is weak and the chaos is strong enough. This limiting case 
should be examined first, and has not been fuUy investigated in previous studies. 

This paper is organized as foUows: In Sec. we introduce a quantum kicked top and its classical counterpart. 
Although this system is well-known in the literature, we again mention them for this paper to be self-contained. In 
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Sec. [II, we introduce the coupled kicked tops (CKTs). We also introduce the von Neumann and linear entropies of 
the subsystem as measures of entanglement. We numerically find that, for both von Neumann and linear entropies, 
the time variation is nearly a linear function of time when the nonlinearity parameter k is large enough. We also find 
that, when the coupling is weak, the production ratc of the linear entropy is nearly proportional to where e is the 
interaction strength between two tops. This result implics that a perturbative treatment is possible. In Sec. IV, we 
derive a perturbative expression for the linear entropy of CKTs. It consists of a correlation function for a single kicked 
top which decays rapidly when the kicked top is chaotic. We also compare numerical results with the perturbative 
expression, and show that the agreement is good up to the long time where the entanglement production rate can be 
defined. In Sec. 0, using the formula in the previous section, we show that the increment of the strength of chaos 
does not enhance the entanglement production rate in the strongly chaotic region. We also numerically confirm that 
this phcnomenon actually happens for CKTs |13 . The relationship between our resul t an d previous results for the 



weakly chaotic region is discussed in Sec. VI. Finally, we summarize this paper in Sec. Vil. 



II. QUANTUM AND CLASSICAL KICKED TOP 

A kicked top is described by the foUowing Hamiltonian: 

n 

where (Jx, Jy, Jz) are angular momentum operators that satisfy [Jx^ Jy] = iJz etc, j is their magnitude, which is a 
conserved quantity, and k is the nonlinear parameter. Since we take 1/j plays an effective Planck's constant. 

The nonlinear parameter k changes the qualitative nature of the classical dynamics (see Eqs. (^) below): On one hand, 
the classical kicked top exhibits chaotic behavior, i.e. the phase space of the classical top is dominated by chaotic 
seas, when k > 3. On the other hand, the classical phase space is dominated by tori when k < 2.5 |Q. A Floquet 
operator (i.e. a one-step time-evolution operator) corresponding to the Hamiltonian (|^) is 

U = e^p{-ikJ^/j) e^p{~ÍTrJy/2). (2) 

In the numerical computations, we employ the |jm)-representation, where J^|jto) = j{j + í)\j'm) and Jz\jm) = 
m\jm). The |jm)-representation of U (Eq. (H)) is (note that j is a conserved quantity), 



where d^'m(/^) Wigner rotational matrix |15|, ll( 



^Ü) ^ -^y-m+m' + m)\{j -my.jj + ni'y.{j -m'y. 



^ (j + m- iy.ll{j - l - m'y.{l - m + m'y. 



X cos(2j'^2'+"-'"')(/3/2) sin(2'-™+™')(/3/2). (4) 

In the numerical evaluation of the Wigner matrix, we employ its Jacobi polynomial expression jï^ . 
The corresponding classical dynamics is described by the foUowing mapping 

x' = z cos{kx) + y sm{kx) , (5a) 
y' = — z siii{kx) + y cos{kx) , (5b) 
z' — —X, (5c) 

where x = Jx/j, y — Jy/j, sjid z — Jz/ j- Note that a;^ + + is a conserved quantity, since j is conserved. We 
hence employ polar coordinates 9 — cos^^ z and (p = tan~^(ï//a;) to concisely describe both classical and quantum 
dynamics. 

In studying the correspondence between quantum and classical dynamics, it is useful to employ a spin-coherent 
State \e,<f) (it) 



{jm\e,(j)) = ————- J— -— (6) 

(1 + |7|^)J V (j +m)!(j -m)! 
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where 7 = e*'^ tan(6'/2). We accordingly employ the Husimi representation of a state vector \^): 

Q{0,cj,)^\{9,m\^- (7) 

To study the classical counterpart of the spm-coherent state \9o,<j)o), we employ a "Gaussian" distribution function 
on the classical phase space (6*, 

foie, </)) ~ e-^o-eo)V2a--i^-M'/2<^' (g) 

where the fluctuation a = is determined by the Husimi function of j^Oj 0o)- 

Figures ^ and ^ show the classical dynamics of the distribution function and quantum dynamics of the Husimi 
function for the kicked top in a semiclassical regime (j = 80), respectively. As is mentioned above, at the initial stage 
of the dynamics, the quantum-classical correspondence holds well. That is, both quantum and classical distribution 
functions bchave very similarly. Such a precise correspondence between the distribution functions is lost as time 
elapses, due to quantum interference. However, even in a much longer time period, the variances of z{t) — Jz{t)/i, 

aUt) = {z{tf),,-{z{t))l, (9) 
^quW - {z{tf)^n-{z{t))l^, (10) 

have a good quantum-classical correspondence. See Fig. ^. Here (. . .)ci and (. . .)qu are the expectation vàlues for the 
classical and the quantum systems, respectively. An explanation in terms of the classical phase-space dynamics is as 
foUows: In the classically regular region (fc = 0.5,1.0), trajectòries are trapped by tori, and the variances exhibits 
periòdic modulations. According to the sizes (along the Jz direction) of the trapping tori, the variances takes various 
vàlues. Furthermore, the variances exhibit recurrence phenomena within a rather short time period (not shown here). 
On the other hand, in the classically chaotic region (fc — 3.0), both (t^[(í) and (Tqy(í) increase rather quickly and 
saturate to the value a"^^^ = 1/3 which is estimated by assuming a uniform distribution. In other words, the phase- 
space distribution functions spread uniformly all over the whole phase space, which is bounded in the case of the 
top. 



III. NUMERICAL EXPERIMENT OF QUANTUM KICKED TOPS 



Here we return to our motivation: How does classical chaotic dynamics affect the entanglement production of the 
corresponding quantum system? To investigate this issue, we employ coupled kicked tops (CKTs), which is introduced 
by Miller and Sarkar ||ï^, as a target system of numerical experiments. The CKTs are described by the following 
Hamiltonian: 

H{t)=Hi{t)+H2Ít)+H,{t) (11) 

where 



H,{t) = -Jz,Jz,Y.5{t-n), (14) 



with [Jxi, Jy-m] — iJzi^im etc. (/, TTi — 1,2). Hcrc ki is the nonlinear parameter of the l-ih top, and e is the strength 
of the coupling between these two tops. Corresponding to the Hamiltonian, Eq. ([ï^), we employ a Floquet operator 
(a one-step time-evolution operator) 

U=U,UiU2 (15) 
where Ui = e-'^^-^U'^^ ^-^^Jyji ^ ^ ^-^k^JlJ2j ^-^^J^J2 ^ ^-^^J.,J.J2_ 

Since we will consider only the case where the density operator of the total system at time í, p(í), describes a pure 
state, we employ entropies of subsystems as measures of quantum entanglement JlSl . More precisely, we employ von 
Neumann and linear entropies of the first top: 

5vn(í) - -Tri{p(i)(í)lnp(iHí)}, (16) 
5n„(í) = l-Tri{pW(í)2}, (17) 
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FIG. 1: Classical dynamics of the distribution function with k = 3.0. (a), (b), (p) corresponds to t = 0, 1, . . . , 15, 
respectively. The center of the initial "Gaussian" distribution function (width cr = 0.1) is located at (0o,0o) = (0.89,0.63). 



where p^^\t) = Tr2{p(í)} is the reduced density operator for the first top. Note that the von Neumann entropy (or 
the linear entropy) of the second top takes the same value as that of the first top, when the whole system is in a pure 
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FIG. 2: Quantum dynamics of the Husimi function with j = 80 and k = 3.0. (a), (b), . . . , (p) corresponds to í = 0, 1, . . . , 15, 
respectively. The center of the initial spin-coherent state is located at (^o, 0o) = (0.89, 0.63). 

state. To calculate Svf^{t) and S'nií (í) numerically, we use the eigenvalues A,(í) of p^^\t) as 

Svn(í) = -^A.(í)lnA.(í), (18) 

i 

S,i„(í) = 1-^{A,(Í)}2. (19) 

i 

We numerically examine the productions of quantum entanglement, using separable states as initial states. In 
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FIG. 3: Time evolution of the variances for (a) the classical and (b) quantum kicked tops. The dotted, dashed, and sòlid lines 
correspond to fc = 0.5, 1.0, and 3.0, respectively. The initial states of (a) and (b) are the same as in Figs. ^ and ^ respectively. 
Since the regular evolutions (fc = 0.5, 1.0) are trapped by tori whose sizes are large in the Jz direction, the variances takes 
larger value than that for the chaotic case (fc = 3.0). 

particular, we focus on the the system parameter dependence, i.e. fci, ^2 and e dependence, of the entanglement 
productions (measured by the entropies of the subsystems) . For simplicity, we only show the cases where k = ki = fc2 . 
As an initial state, we employ the foUowing pure and separable state 

p(0) = |vl/(0))(*(0)|, (20) 

where |í'(0)) = \9i, 0i) (g) |02, ^2) is a direct product of spin coherent states. In studying the chaotic region, the center 
of the initial spin coherent state is located in the chaotic sea. Even such numerical experiments with a restricted class 
of initial states provides an insight about the typical behavior of chaotic CKTs, when the fraction of tori is small in 
phase space for the corresponding classical system (in the case of CKTs, k > 3.0). 

Figure ^ shows the time evolutions of the entropies with various vàlues of fc. The entropies stick to nearly zero 
until the "rising time", and then they increase nearly linearly as a function of time for chaotic cases (see Figs. ^ 
(b),(c),(d)). Though not shown here, they saturate to finite vàlues after the long time evolution, due to the finiteness 
of the dimension of the Hilbert space. In the foUowing, we focus on the "intermcdiate" region where the entropies 
increase monotonically as a function of time. Our extensive numerical experiments conclude that, in the intermcdiate 
region, the linear entropy as well as the von Neumann entropy increases nearly linearly as a function of time, when 
the chaos is strong enough and the coupling is weak enough. 

Figure ^ shows e dependence of the linear and von Neumann entropies at T = 128. When the coupling is weak 
(i.e., e < 10^'^), the linear entropy obeys a perturbative behavior S'iin(r) cx (see Fig. ||(a)). At the same time, the 
time step T — 128 belong to the region where the linear entropy increases linearly in time. In contrast to this, when 
e is much larger than 10"'^, the entropy does not belong to the perturbative region, and saturates to a finite value, 
which is determined by the finite size of the Hilbert space of CTKs. These observations suggests us to analyze the 
behavior of the linear entropy using a perturbation treatment for the interaction strength e. This is the subject of the 
next section. On the other hand, we confirmed that S'vnÍ^) oc e"'^'* when e is small enough (Fig. || (b)). It seems that 
an usual perturbative treatment is dif&cult to explain the exponent 1.8, so we will concentrate on the linear entropy 
below. 



IV. PERTURBATIVE EXPRESSION FOR THE LINEAR ENTROPY 

A. Perturbation treatment 

We evaluate 5'iin(í), Eq. (p7|), the linear entropy of the first top, by using the time- dependent perturbation theory 
with a small parameter e. First, we introduce the interaction pictures, of the density matrix p{t) = C/q (í)p(í)/7o(í), 
and, of the operator A{t) = UQ{t)AUo{t), where C/o(í) = (Í7i<8ií72)*- That is, A{t) corresponds to the "free" evolution of 
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/ t 

FIG. 4: Time evolutions of the linear (dashed line) and von Neumann (soIid line) entropies. Parameters are e — 10"'' (in an weak 
coupling region) and j = 80 (in a semiclassical region). The initial state is Eq. (BOI) with 9i — 02 — 0.89 and 4>i = (t>2 = 0.63. 
The center of the initial state is located in the chaotic seas for the cases of (b), (c), and (d). 
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FIG. 5: e dependence of (a) the linear and (b) von Neumann entropies at T = 128 with k — 3.0. The dashed lines denote 10^ e'^ 
and 10''e^ ** for (a) and (b), respectively. The initial state is the same as in Fig. 0. 
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the operator A. Accordingly, the time evolution of p{t) is described by the unitary mapping p{t) = Ü^{t)p{t— l)í7j'(í), 
where the expansion of Uc{t) by small e takes the following form 



2 



U,{t)^l--V{t)-—V\t)+0{e^) (21) 
with V{t) — J zi{i) J z2{i) ■ Hcncc, the unitary mapping of p{t) becomes 



,2 



m - p{t - 1) + -[p{t - l),V{t)] - —^[[p{t - 1), í/(í)], V(t)] + 0{e'). (22) 



By induction, we have 

pit) = p(o) + -^[p(o),no] 



t 

te 



^ 1=1 

2 * 2 * 

4 E E [[^(0)' ^(01. ^M] - L· E[[^(0)' ^(01. V^(0] + 0{e^). (23) 

/=2 m=l 1=1 



By tracing out the second system, we have 

í 

pW(í) = p«(o) + ^^[pW(o),j,i(0](J.2(0>2 



àE{[^-i(Op^'^(o)'·^^iW] + [·^^i(0,p^'Ho)j,i(0]}(ifi(0>. 



1=1 

g2 t í~l 
/=2 m=l 



1=2 m=l 

+ [J.i(m)p(i)(0), J.i(0](Í.2(0^.2(m))2} + O(e') (24) 
where (A)2 = Traíp^^) (0)^} is an average for subsystem 2. Finally, we obtain a second order perturbation formula of 

t t 

5L^(í) = 5oEE^(^'"^) (25) 

1=1 m=l 

where Sq = 2e^j^ and D{l,m) is a correlation function of the uncoupled system. Since the interaction Hamiltonian 
Hg{t), Eq. (^4|), is in a bihnear form, D{l,m) is decomposed into a product of correlation functions of uncoupled 
subsystems 

D{l,m) Ci{l,m)C2{l,m) (26) 

where 

Cdl.m) = (z;(Of.(m)), - (f,(0)(i.(m)), (27) 



and Ziiï) = Jzi{l)/j (i = 1, 2). In the perturbation formula, Eq. (p5|), S'ò is a rather trivial factor implying S{'^J{t) oc 
as is observed in Fig. |^(a). The nature of the dynamics for the tops is reflected in D{l,m). 

Let us remark important points of our perturbation formula: (i) In common with the exact case, S^^J{t) has a 
symmetric form about the exchange of the first and the second tops. That is, our perturbative treatment preserves 
this symmetry, although we start from a perturbative treatment of the linear entropy of the first top. (ii) Our 
formula has a similarity with those in phenomenological descriptions of linear irreversible processes ||20|] , in the sense 
that these theories use time correlation functions to describe relaxation phcnomena. This is useful both for making 
phenomenological arguments and for establishing a link between a phenomenological theory and a microscopic theory 
(cf . the linear response theory of nonequilibrium statistical mechanics [ pO[ ) ; (iii) Since our approach does not take into 
account the effect of the recurrence, the formula ( p5| ) would have qualitatively different applicability to the classically 
regular and chaotic systems. For classically regular systems, our theory would break down in relatively short time 
period, due to the smallness of the period of the recurrence. On the other hand, for chaotic systems, we numerically 
confirmed that our theory works for a rather long time period. 
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FIG. 6: Estimations of linear entropies by the perturbative formula, Eq. (ESI). It is compared with the exact numerical result 
(in dashed lines). The parameters are (a) e = 10~* (perturbative regionjwith k = 1 (regular) (b) e = 10~* (perturbative 
region) with k — S (chaotic), (c) e = 10"'^ (intermediate region) with k — 1 (regular), (d) e = 10~^ (intermediate region) with 
k = 3 (chaotic). 



B. Comparison with numerical results 

We numerically examinc our formula, Eq. (p5|). In Fig. ^ we plot both S'nií (í) and S^iJ{t) for the intermediate 
coupling and weak coupling cases with regular and chaotic conditions. The initial state is the direct product of the 
spin-coherent state as before. As shown in Fig. |^, e = 10~^ is the perturbative region, so the agreement between 
Sun{t) and Sf^J{t) is very good for different k's up to í ~ 100 (Figs. ^ (a) and (b)). Note that our perturbative 
expression works for such a long time to reproduce the linear increment of the entanglement productions in time. 
Such a correspondence degrades as e gets larger, of course, as shown in Figs. || (c) and (d). However, as far as 
concerning the chaotic case k = 3.0, our expression describes the entanglement production, at least, qualitatively. 



V. DYNAMICAL ASPECTS OF ENTANGLEMENT 

A. Harder chaos does not mean larger entanglement production 

In this section, the perturbative formula, Eq. (p5|), is employed to answer the foUowing question: how does the 
strength of chaos influence on the entanglement production rate in the strongly chaotic regions where the influences 
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from tori are negligible. 

We examine Ci{l,m), Eq. (p7|), which describes the fluctuation of Zi — Jzi/j- Since the tops are strongly chaotic, 
we impose several phenomenological assumptions on Ci{l, m). Since the phase space of the kicked top is bounded, the 
distribution function in the phase space becomes quickly uniform in the strongly chaotic region. Hence we assume 
— "'satJ where we ignore a short transient before the distribution function becomes uniform (see Fig. |3|). The 
magnitude of the fluctuation cr^g^^ = 1/3 is determined by the assumption that the distribution function is uniform 
on sphere {0, (j)). The boundedness of the phase space allows us to employ another assumption that the relaxation of 
Ci{l,m) [l ^ m) is exponential with an exponent 7^ [ p^ . Furthermore, it is natural to assume that the exponent 7^ 
becomes larger as the positive Lyapunov exponent of the corresponding classical system becomes larger. The simplest 
function that satisfies the assumptions above is 

C.(/,m)^aLe-^'l'-™l. (28) 

Hence D{l,m), Eq. (|2^), becomes 

i:)(/,m) ~ Doe"^!'""! (29) 
where Dq = a^^^ and 7 = 71 + 72 ■ Accordingly, Eq. ( p5| ) provides the foUowing evaluation of the linear entropy 

coth(7/2)í 1 ~ e 



^lin (i) — SqDq 



sinh 7 — 1 



(30) 



Whcn the relaxation time of D(l, m) is much shorter than the time scale of the stationary entanglement production 
region, we have an entanglement production rate T 



dt 



~roCoth(7/2), (31) 

í>l/7 

where Fq = SqDq. From this relation, it is shown that F decreases as 7 becomes larger, i.e. the chaos of the 
corresponding classical system becomes stronger. That is, the increment of the strength of chaos does not enhance 
the production rate of entanglement. Furthermore, in the limit 7 — > 00, F quickly saturates to a finite value, Fq. 

At first glance, our prediction seem to be counter-intuitive. Hence we provide an explanation of the prediction to 
summarize this subsection: The entanglement productions are induced by the fluctuation of the intcraction Hamil- 
tonian H^{t), in the interaction picture. Since the time dependence of H^{t) looks like very "random" in classically 
chaotic Systems, the contribution from H^(t) to the linear entropy is reduced due to dynamical averagin n. We note 
that this mechanism is similar to that of the so called motional narrowing in spin relaxation phenomena po| , ^ï]| . 

B. Properties of the correlation function and linear entropy production rate — Saturation of the 

entanglement production 

We test the prediction of the phenomenological argument above with numerical experiments. First, we examine 
the correlation function D, Eq. ( p6| ) (see Fig. We confirmed the assumption, Eq. (p9|), for D when the classical 
counterpart is chaotic (fc > 3.0): The correlation function decays very quickly, althogh it is difhcult to detemine the 
exponent 7 directly from the numerical evaluation of D (see another estimation of 7 below). On the other hand, the 
value of D{t,t) is almost independent with k, and approximately equal to Dç, — a"^^^ = 1/9. 



Second, we examine the nonlinear parameter fc-dependence of the entanglement production rate F, Eq. (31). For 
each initial states, whose center of the spin-coherent states is placed in the chaotic sea, we obtain F using the least 
square fitting for the time region from í = 20 to 100, where í-linearity holds (Fig. |^). Hence we confirm that the 
increment of the strength of chaos does not enhance entanglement production rate in the perturbative regime, where 
e is small enough. It is also confirmed that the entropy production rate F saturates to Fq for large k, which is also 
consistent with Eq. (^). At the same time, we numerically find that Fvn, the entanglement production rate measured 
by von Neumann entropy also exhibits a saturation in the large k limit. Although we do not have any analytical 
theory for Fvn, we cxpect that the saturation of F^n is explained by the similar explanation as that of the linear 



entropy F (see Sec. VA ). 

In Fig. H we plot fc-dependences of two quantities: One is Asum = Ai + A2, where \i is the short-time (up to í = 100) 
and phase-space averaged Lyapunov exponent for the initial classical distribution of i-th subsystem. We note that 
Asum = 2 Al since fci = ^2- The other is the decay rate 7 of the correlation function D, Eq. (p6|). We estimate 7 from 
the phenomenological estimation of F, Eq. (^I]), i.e. 

T/V 
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FIG. 7: Plot of D^{t,t~ t)/Do (= ^D{t,t - r)/Do) as a function of r with different í's (í = 40,50,60,70) and ifc's. (a), (b), 
(c), (d) corresponds to = 1, 3, 5, 7, respectively. The initial state is Eq. ( po| ) with 6i = 62 ~ 0.89 and 4>i ~ 4'2 ~ 0.63. 

instead of the direct estimation from thc assumption, Eq. (p9[). Figure ^suggests 7 ~ Asum- This shows an evidence 
that the decay rates of the correlations of tops are determined by the positive Lyapunov exponents of the classical 
counterparts. Thus our numerical experiments confirm the estimation, Eq. (pT|), and it is concluded that the entan- 
glement production rate is not increased by the increment of the strength of chaos in the strongly chaotic region, and 
saturate to a finite value in the strong chaos hmit. 

Finally, we point out that it is natural to generahze our study on CKTs to any strongly chaotic system with bounded 
phase space. We note that we have already confirmed this for coupled kicked rotors [|[ |ï^, with periòdic boundary 
conditions of both position and momentum coordinates. 



C. An extension to flow systems 



In this subsection, we extend the above argument to flow systems with continuous time. Consider the casc where 
flow systems are weakly coupled. When thc initial state is a pure product state (as is the case above), the linear 
entropy produced in the composite system is 



So í dT [ dr'Diry), 
Jo Jo 



(33) 
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FIG. 8: (a) k dependence of the normalized linear entropy production rate for various initial conditions located in the chaotic 
sea. Note that F saturates to To as k increases. (b) k dependence of the normalized von Neumann entropy production rate Fvn 
for various initial conditions (the same as those of (a)) located in the chaotic sea. Note that Fvn is scaled by Fo = 2e^'^j^Do 
instead of Fo = 2e^fDo (cf. Fig. |). 




FIG. 9: k dependence of the correlation decay rate 7, which is determined from Eq. (p2|), for various initial conditions located in 
the chaotic sea (scattered symbols). The sum of the Lyapunov exponent Xsum is also shown as a bold curve, which correspond 
to a single initial condition. We note that Asum depends on the initial condition only very weakly. 



where D{t, t') is a correlation function determined by the form of the interaction Hamiltonian. When the subsystems 
are strongly chaotic with bounded phase space, we assume again 

D{t,t') ^ Doe-^\^-^\ (34) 

Substituting this into Eq. we have 

r^?Míl(i_e-7*)_^MíI (í»l/7). (35) 
7 7 

Hence, if the subsystem is strongly chaotic (i.e. 7 — > cxd), the entanglement production rate becomes zero. That 
is, strong chaos completely suppresses entanglement production! Although we have not numerically confirmed this 
suppression yet, the similar suppressions of quantum relaxations in str ongl y chaotic systems have been observed by 



Prosen and Znidaric |2a, M. We will discuss this point further in Sec. VII 
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VI. DISCUSSION ON A WEAKLY CHAOTIC REGION 



In the previous section, we investigated the strongly chaotic region of thc CKTs and concluded that the increment 
of the strength of chaos does not enhance entanglement productions, which is measured by the hnear entropy of the 
subsystem, with the help of the perturbative formula, Eq. (^5|). With this in mind, we discuss the recent study by 
MiUer and Sarkar , who investigated the weakly chaotic region (where chaotic seas and tori coexist) of the CKTs, 
and claimed that the increment of the strength of chaos enhances entanglement. More precisely, they numerically 
found that the entanglement production rate, which is measured by the von Neumann entropy, linearly depends on 
the sum of positive (finite-time) Lyapunov exponents of the corresponding classical CKTs, without any theoretical 
justification. As is well known, it is much harder to develop a theory of weakly chaotic systems (in other words, mixed 
phase space systems) than strongly chaotic systems. This is actually the case with the numerical result of Miller and 
Sarkar. To accommodate these two qualitatively different results, we employ our perturbative formula, Eq. (|2^), in 
the analysis for the weakly chaotic region (fc = 3.0, Fig. [l·Ol(a)). 

In order to justify the application of our formula, Eq. (|2q), we confirm that the entanglement productions measured 
by the linear entropy, instead of the von Neumann entropy, reproduce Miller and Sarkar's fitting. See Fig. [ÏO|(b). 
Furthermore, we numerically examined that the perturbative evaluation of the entanglement production rate F, 
Eq. (|2^), is applicablc to the weakly chaotic regions, in particular the case above. Hence, in the foUowing, we 
reexamine the inputs of the formula, Eq. (p5[), which is the correlation functions of the uncoupled systems. 

We focus on our assumption, Eq. (p9|), for the correlation functions Z3(Z,to), Eq. (^), which is derived from 
Eq. ( psi) for strongly chaotic regions: (i) Due to the absence of tori in thc corresponding classical system, the 
fluctuation of Zi takes a saturated value which agrees with that of the uniform distribution in the classical phase space 
(i.e., a{l,l) 

— '''sat); (ii) -Due to the strongly chaotic dynamics, the correlation functions decays exponentially (i.e. 
Ci{l,m) oc e-T-l'-"!). 

First, we examine our assumption that a2, the fluctuation of Z2 — Jz/ ii is independent of 62. In the weakly chaotic 
region, this assumption breaks down due to the confincment of phase-space dynamics by tori. Actually, (T2 becomes 
smaller as the "overlapping" between the initial state and tori become larger (see Fig. ^(a)). By taking account of 
this fact into the assumption Eq. ( p9[ ) on DÜ,m), the decrement of (T2 in tori provides a crudest explanation of the 
6*2 dependence of F (denoted by □ in Fig. |l^(b)). That is, the decrement of the fluctuation of Z2 (denoted by □ in 
Fig. ^ (a)) due to the influence from tori inhibits the entanglement production. See the line with □ in Fig. ^ (b). 
However, the improved estimation denoted by in Fig. |ïï| (b) still exhibits a qualitative discrepancy. 

Second, to overcome this discrepancy, we improve the assumption for D{l,m) as 

D{l,m) = a2cr2e-^l'-"le'''^('-'"\ (36) 

where we introduce a real-valued parameter to. This characterizes oscillations due to the regular motion of the second 
top. The resultant oscillation of D{1, m) tends to reduce the value of F in Eq. (p5| ) (cf. Eq. (^) : 

-p (o-l/o-sat)^ (o·2/a·sat)^ -p 

^ = 1 I r ■ / loM ■ \.í /0M2 ^ To COth 7/2). 37) 
1 + {sm(ti;/2)/ smh(7/2)}^ 



We determine the value of lü in Eq. (36) from the Fourier transformation of D{l,m) (see Fig. |l·lKa,)). We depict 
the estimation Eq. ( ^7| ) (denoted by A) also in Fig. |ïï| (b). We conclude that the assumption Eq. ( |3q ) provides a 
satisfactory improvement of the evaluation of F for the weakly chaotic region that Miller and Sarkar investigated. 

From our argument, it is seen that the contribution from tori also play a role for the determination of the entangle- 
ment production rate via 172 and lü. Thus it is suggested that the linear dependence of F with the sum of the positive 
Lyapunov exponents of the corresponding classical system is not intrinsic for the weakly coupling region. 



VII. SUMMARY AND OUTLOOK 



We have studied how the strength of chaos affects the production rate of quantum entanglement of the coupled kicked 
tops (CKTs). When the coupling constant e is small enough, the entanglement productions obey the perturbative 
formula, Eq. (|2^) . When the classical counterpart exhibits chaotic behavior, there appears a "stationary" entanglement 
production regime where the entanglement production rate is well-defined. In the strongly chaotic limit, where the 
correlation functions of the uncoupled tops decay exponentially fast, the perturbative formula, Eq. (p5[), predicts that 
the entanglement production rate saturates to a certain value. Our numerical experiment confirmed this prediction. 
This is a unexpected result since the previous works show that the chaotic dynamics promotes a larger amount of 
quantum entanglement compared with the regular dynamics 0, |[ |[ |Ï0) • 



14 




FIG. 10: (a) A magnification of Poincare section of the kicked top with k — 3.0. (b) The correlation between Asum and F in 
a weakly chaotic region. In particular, sòlid circles correspond to those in (a). A linear fitting is depicted by a dashed line 
(r/Fo — aAsuni + b with a — 2.2 and b = —0.84). Parameters are e — 10~^ (perturbative region) and j — 80. The center of the 
initial state of the first top is located at (çii, ^i) = (0.63, 0.89), and those of the second top are depicted by sòlid circles in (a). 




FIG. 11: (a) 02 dependences of (0'2/o'Buin)'^ (□), Aaum (O)) lo (A). (b) The 62 dependence of F/Fo (□). The perturbative 
estimations with the assumptions (09) and lam are indicated by Q s^iid ^) respectively. Parameters and initial conditions are 
the same as in Fig. (b). At the same time, we assume (cri/íTsat)^ = 1 for the perturbative estimations. 



Our perturbative argument of the strongly chaotic region depends only on the two points: (i) the time-correlation 
function of the interaction Hamiltonian decays exponentially; (ii) the phase space distributions of the corresponding 
classical subsystems become quickly uniform before the stationary entanglement production starts. Hence we expect 
that our result also holds for a wide variety of classically chaotic systems. 

At the same time, we reexamined the weakly chaotic region which is recently investigated by Miller and Sarkar, who 
showed numerically that the entanglement production rate linearly depends on the sum of the positive linear stabiUty 
exponents Our perturbative approach provides a theoretical way to explain their result: The entanglement 

production rate is controUed by the combination of the decay rate, the magnitude and the oscillation frequency of 
the time-correlation function of the interaction Hamiltonian. It is hard to believe that these factors are generally 
determined only by the Lyapunov exponents. Rather, it is natural to expect that the behavior of the correlation 
function is strongly influenced by the existence of tori. 

We point out that the investigation of dynamical production of quantum entanglement has relevance with that of 
quantum fidelity (measured by an overlapping integral of two states that are evolved by slightly different Hamiltonians) 
p^ , p5| . The decay of fidelity and the production of entanglement correspond to the quantum relaxations against 
static and dynamic disturbances, respectively. When the disturbance is small enough, the perturbative approach will 
describe the leading ("linear") response. We showed that this is the case for the dynamical production of quantum 
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entanglement. On the other hand, conceriiing the evaluations of quantum fidelity, Prosen et al. reported the success 
of a perturbative approach ||25|| . Both works predict that the strong chaos suppresses the quantum relaxations in 
flow Systems. Furthermore, Prosen et al. reported that their theoretical prediction on the static disturbances agrees 
with their numerical experiments ||2^ . The recent studies of quantum fidehty for nonperturbative regimes [^6| will be 
apphcable to the studies on dynamical productions of quantum entanglement. We beheve that such an effort will be 
fruitful to investigate "quantum chaos" in many degrees of freedom systems (see, e.g., Refs. ^ï\). 

Finally, we point out a possible application of our work to the studies of realistic systems. In the investigations of 
Chemical systems with large degrees of freedom including biological systems [^9| , it is important to estimate the 
entanglement (decoherence) rate. Most studies on this problem rely on the approaches using the master equations 
or the influence functional technique. However, these approaches have a serious difhculty in practical applications 
to Chemical reaction dynamics, since the time-scale separation of the two constituents ("the system" and "the en- 
vironment") in the whole system often breaks down. In contrast to this, our approach only assumes the weakness 
of the coupling between the subsystems, which dynamically causes entanglement in the whole system. Hence it has 
an ability to cope with the breakdown of the the time-scale separation. We expect that our approach will provide a 
useful tool to investigate chemical reaction dynamics. 
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